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Abstract 



We study the thermodynamic properties of spin chains of Haldane-Shastry type associated with 
the Ajv-i root system in the presence of a uniform external magnetic field. To this end, we exactly 
compute the partition function of these models for an arbitrary finite number of spins. We then 
show that these chains are equivalent to a suitable inhomogeneous classical Ising model in a 
spatially dependent magnetic field, generalizing the results of Basu-Mallick et al. for the zero 
magnetic field case. Using the standard transfer matrix approach, we are able to compute in 
closed form the free energy per site in the thermodynamic limit. We perform a detailed analysis 
of the chains' thermodynamics in a unified way, with special emphasis on the zero field and zero 
temperature limits. Finally, we provide a novel interpretation of the thermodynamic quantities 
of spin chains of Haldane-Shastry type as weighted averages of the analogous quantities over an 
ensemble of classical Ising models. 

Keywords: Spin chains of Haldane-Shastry type, vertex models, transfer matrix method, 
thermodynamic limit 



1. Introduction 

In this paper we study a class of su(m) spin chains whose Hamiltonian can be collectively 
written as 



where the B a \ are real constants and the interactions J t j are chosen as described below (see 
Eqs. (1.7)). In the previous formula, the operators S t j act on a state 




(l.D 



\s u ...,s N ), 



of the canonical spin basis by permuting the i-th and j-th spins: 



Sij\su. ..,s h ...,Sj,..., s N ) = \su- ,J,-, ... , s N ) . 
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T? with the normalization \x(T?T^) - 26 a/3 (where the subindex k labels the chain sites) as [ 1 1 



The permutation operators S y can be expressed in terms of the (Hermitian) su(ra) spin operators 

m 2 -l 

2 



1 1 " r ~' 

+ (1-2) 



<r=\ 



Let S k denote the operator whose action on the canonical spin basis is given by 

Sl\s u ...,s N ) = (6 a Sk -Sl)\s 1 ,...,s N ), a=\,...,m-\, (1.3) 

so that the operators LS? span the standard Cartan subalgebra of su(ra). The operators S a are 
then defined by 

N 

S a = YS?, a= 1 m- 1. (1.4) 

(=i 

Thus the second sum in Eq. ( |1.1[) c an be interpreted as arising from the interaction with a uniform 
external su(m) "magnetic" fieldMof strengths (B\, . . . ,B m _i). Note that in the case m = 2 (i.e., 
for spin 1/2) we can take {7^, T^, T k ) = {cr x k , cr y k , cr^} and S| = o\, where <T v k is a Pauli matrix 
acting on the A:-th spin's Hilbert space. Hence Eq. \1.2) adopts the more familiar form 

S ij =\{l+o- i -<Tj), (1.5) 

and the Hamiltonian ( |1.1[ ) reduces to 

<H = J] ^(i-^.^-B^., (1.6) 

with B = B\. In particular, the last term represents the interaction with a uniform magnetic field 
parallel to the z axis with strength (proportional to) B. 

The three models we shall study are defined by the following choice of the interaction 
strengths J if 

• The Haldane-Shastry (HS) chain |2|3): 



/ _kn 
2 sin 2 (£ -£,■)' 



The Polychronakos-Frahm (PF) chain |U|3): 



■/;,= — (l-7b) 

J (£-ft) 2 



'The most general su(m) magnetic field term is of the form T = where 

m 2 -l 

T, = ]T b " T i 

a=l 

is a traceless Hermitian matrix acting on the internal space of the i'-th spin. By performing a rotation in this internal space 
we can diagonalize the matrix 7j, effectively replacing it by a traceless diagonal matrix. The latter matrix can in turn be 
expressed in the form EJjTj B a S° , which yields the last term in Eq. \\.\) . 
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where & is the &-th root of the Hermite polynomial of degree N. 
• The Frahm-Inozemtsev (FI) chain 1 6 1 : 



J i} = ~ 7 , (1.7c) 

2sinh 2 (£-^) 

where e 2 ^' is the k-th root of the generalized Laguerre polynomial L a N x with a > 0. In all 
three cases, J is a real constant whose sign determines the model's ferromagnetic (J > 0) or 
antiferromagnetic (J < 0) character. Note that, while the sites of the HS chain are equispacecQ 
this is not the case for the PF or FI chains. 

We shall denote by 

<H = J] (1.8) 

l<!<7<iV 

the Hamiltonian of the chains in the absence of a magnetic field. Following standard ter- 
minology, we shall collectively refer to the chains ( |1.7[ )-( fTT8] l as spin chains of Haldane-Shastry 
type. They are all associated with the root system A#_i, in the sense that the interactions Jjj 
depend only on the differences of the site coordinates Although several generalizations of 
these chains to the BCn and root systems have been considered in the literature I7 HTTI . in 
this paper we shall restrict ourselves to the above A^v_i -type models. 

Spin chains of HS type are the simplest models in condensed matter physics exhibiting frac- 
tional statistics 1 12]. Historically, the HS chain ( |1.7a| i-( [TT8"j ) was introduced as a simplified ver- 
sion of the one-dimensional Hubbard model with long-range hopping, from which it can be 
obtained in the half-filling regime when the on-site interaction tends to infinity fOl . Soon after 
its introduction, it was shown that this chain is completely integrable, in the sense that one can 
explicitly construct N -I mutually commuting integrals of motion lfl4l[T5l . As first observed by 
Polychronakos (4), these integrals of motion can be obtained from those of the dynamical spin 
Sutherland model [16| by means of the so-called "freezing trick". In fact, the freezing trick can 
also be applied to derive the PF and FI chains from the Calogero [ 17 1 and Inozemtsev IfTSl spin 
dynamical models. In particular, these two chains are also completely integrable. Apart from 
their integrable character, spin chains of HS type appear in many areas of current interest in both 
physics and mathematics, such as quantum chaos |fT9"ll2"0l , supersymmetry II2TII221 . conformal 
field theory G3TI251 . the AdS-CFT correspondence l26ll27i l. one-dimensional anyons 11281 and 
Yangian quantum groups [15 29 -31]. 

The partition functions of all three chains of HS type in the absence of a magnetic field, which 
have been computed in closed form using again the freezing trick 0214341 . can be expressed in a 
unified way by the formula 

Zo(q) = 2 fl (TV) ■ <r mK,) R (1 " " (L9) 

ke<P N i=l i=\ 



2 In fact, if 2^ is interpreted as an angular coordinate, then the HS chain describes an array of spins equispaced on 
the unit circle, with long-range pairwise interactions inversely proportional to the square of the chord distance between 
the spins. 
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Here q = e 1 / ( - k B T \ ks (k\, . . . ,k r ) is an element of the set P N of partitions of N with order taken 
into account, and the numbers K h K'j e (1, . . . , N — 1} in Eq. ( |1.9| l are positive integers defined by 

£ 

Ki = Yj k p tK[,...,K' N _ r } = {l,...,N-l}\{K u ...,K r - 1 }. (1.10) 

j=i 

Remarkably, the partition function j£o of the chains ( |1.7| )-( fTT8) l depends on the chain under con- 
sideration only through its dispersion relation 



i(N - i) , for the HS chain 

i, for the PF chain (1.11) 

i(a + i - 1) , for the FI chain. 




£(0 



Using the previous expression for the partition function, Basu-Mallick et al. 1 35 1 showed that 
the spectrum of the spin chains of HS type is given by 

N-l 

E (n) = J^6(n h n M )&(i) , n = (n u . . . , n N ) , (1.12) 

;=i 

where 

j > k 

■ <• v (L13) 

and the quantum numbers «,■ independently take the values 1, . . . ,m. The vectors d(n) e {0, l}^ 1 
with components <5,(n) = <5(n,-,n, + i) are in fact the celebrated motifs introduced by Haldane 
et al. 061 . As pointed out in Ref. 11351 . Eq. ( |1.12| i defines a classical inhomogeneous one- 
dimensional vertex model with bonds n\, . . . ,n^ each taking m possible values, where the con- 
tribution of the k-th vertex is given b^F\j6(nic-\,nk)S(k - 1). We shall show that this connection 
between spin chains of HS type and vertex models still holds in the presence of a nonzero su(m) 
magnetic field. This is indeed the key result on which we shall base our unified approach to the 
analysis of the thermodynamics of all three chains of HS type ( |1.1[ ). 

The study of the thermodynamic properties of spin chains of HS type was initiated by Hal- 
dane himself, who used the spinon description of the spectrum to derive an indirect expression 
for the entropy of the spin 1/2 HS chain [23]. An explicit formula for the free energy of the PF 
chain in the absence of a magnetic field appeared shortly afterwards without proof in Ref. (5). 
In a subsequent publication [6|, Frahm and Inozemtsev obtained an analogous expression for 
the FI chain using the transfer matrix method, and also computed its magnetization for arbitrary 
magnetic field. 

In this paper we have several objectives that we shall now summarize. In the first place, we 
shall compute in closed form the partition function of the HS-type chains (jTTTJ in the presence 
of an arbitrary su(m) magnetic field (B\, . . . ,B,„_i). Secondly, we shall use the expression for 
the partition function to establish the equivalence of the latter chains to a suitable classical in- 
homogeneous vertex model. We shall then take advantage of this connection to determine the 
equilibrium thermodynamics of the spin 1/2 chains of HS type in a unified and systematic way. 



3 By Eq. (1.12) , the first and last vertices do not contribute to the energy. 

4 



Finally, we shall use the previous results to uncover a novel connection between spin chains of 
HS type and the classical Ising model. 

The paper is organized as follows. In Section [2] we compute the chains' partition functions 
in closed form by means of the freezing trick. Following the approach of Ref. Il35ll . in Section[3] 
we construct a generating function for the partition function in terms of complete homogeneous 
symmetric polynomials. We then define a similar function for the corresponding vertex models, 
and prove that both generating functions coincide for arbitrary values of their arguments. This 
shows that the energies of each of the chains (jTTTJ are identical to those of its associated vertex 
model. Exploiting this connection, in Section|4]we use the transfer matrix method to find explicit 
expressions for the thermodynamic functions of the chains (jTTTJ with spin 1/2. In the case of the 
PF chain, we derive an exact formula for the free energy in terms of the dilogarithm function in an 
arbitrary magnetic field, and show that its magnetization can be expressed in terms of elementary 
functions. Section[5]is devoted to a detailed study of the zero magnetic field case. In particular, 
we show that the susceptibility can be expressed for all three chains of HS type in terms of the 
error function, discuss the connection of these chains with two-level systems, and derive low- 
temperature asymptotic expansions the chains' thermodynamic functions. In Section|6]we study 
the zero temperature limit in the presence of an arbitrary external magnetic field in both the 
ferromagnetic and antiferromagnetic regimes. We explicitly show that in both regimes there is a 
phase transition, as was to be expected on general grounds [5|. In Section|7]we present a novel 
interpretation of the thermodynamic quantities of spin chains of HS type as weighted averages 
of the analogous quantities over an ensemble of classical Ising chains. The paper ends with a 
short section where we compare our results with previous related work, and comment on possible 
extensions thereof. 



2. Partition functions 

In this section we shall derive a closed form expression for the partition functions of the HS 
chains (JTTTJ in the presence of a constant su(m) magnetic field. For definiteness, we shall deal 
with the PF chain, whose interactions are defined by Eq. ( |1.7b| i. As we mentioned above, this 
chain is obtained by applying the freezing trick to the su(ra) spin Calogero model in a constant 
magnetic field, whose Hamiltonian we shall take as 

~ m-l 

H = H - — YB a S a , (2.1) 

J a=l 



where a > 1/2 and 



More precisely, let 



;=i ;=i \<<*i<N y 1 J' 



denote the Hamiltonian of the scalar Calogero model 1371 . and define 



1 _ C , m-l 

\^i<j^N v ' J ' a=l 
5 



(2.4) 



so that 

H = H SC + 2a <#(x) , <H = /<#(£) . (2.5) 

It is well known 1 38 1 that the chain sites £ of the PF chain are the coordinates of the unique 
minimum of the scalar potential 

t/(x) = f>f + £ —1-5 (2.6) 
,=1 l^i^JV w '' 

in the principal Weyl chamber of A w _ 1 type jci < x^ < ■ ■ ■ < x N , where the motion of the particles 
of the Calogero model can be restricted due to the singularities of its Hamiltonian. From this fact 
and Eq. ( |2.1[ ) it follows that in the large a limit the eigenfunctions of the spin Calogero model 
become sharply peaked around the sites of the PF chain, so that the su(m) degrees of freedom 
effectively decouple from the dynamical ones. By Eq. ( |2.5| ), in this limit the energies of the 
model ( |2.1[ ) are approximately given by the formula 

E tj *Ef+yej, (2.7) 

where E, and ej denote two arbitrary eigenvalues of H sc and < H, respectively. Although the 
spectra of both the scalar and the spin Calogero models can be easily computed, the previous 
formula cannot be directly used to compute the spectrum of the PF chain. Indeed, it is not 
obvious a priori which energies of H and fP c can be combined to yield an energy of the PF 
chain. However, from Eq. ( |2.7| i it is immediate to obtain the following exact relation between the 
partition functions Z, Z sc and Z. of the Hamiltonians H, H sc and *H: 

Z(T) ^ lim Z(2aT/J) ■ (2.8) 

This formula, first derived by Polychronakos |32|, is the mainstay of the freezing trick method; 
see 1 39 40 1 for a rigorous proof. 

We shall next evaluate both the numerator and the denominator in the previous equation. To 
begin with, the partition function of the scalar Calogero model is easily computed 021 . with the 
result 

JV 

Z sc (2aT/J) = q% ]~[(1 - q u y l , (2.9) 
i=i 

where 

E g =a 2 N(N-l) + aN (2.10) 

is the ground state energy. 

Let us now turn to the partition function of H. Consider, to this end, the spin functions 

fe(x) = p(x)A( flxf s N )) , (2.11) 

i=l 

where l t 6 NU {0}, t = (€ u . . . ,t N ), s = (s u . . . , s N ), 

N 

p(x) = e- 2 / 2 ]~| \ Xi - Xj \ a , r 2 = Y^ x l 

l^i<j%N i=l 
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is the ground state of the scalar Hamiltonian H sc and A is the total symmetrizer with respect to 
particle permutations. The above states are a (non-orthonormal) basis of the Hilbert space of H 
provided that (for instance) ti ^ for all i and s,- > S{+\ whenever t\ = It was shown in 
Ref. [4-1] that Ha acts triangularly on the latter basis, with eigenvalues 

N 

E g +2a\e\, \f\ = Yj€ t . (2.12) 

!=1 

On the other hand, from the identities 

SijSj = SjSij , SijSj - SfSij '< S'jSk = Sl-Sij ' k&hji 

it immediately follows that the operators S a commute with the spin permutation operators S 
and hence with A. By Eq. ( |1.3) we then have 



N 

n 

' i=1 

where the eigenvalue c a (s) is given by 

N 



S> f , s (x) = p(x)A( Uxf S a \s u ..., s N )) = c a (s)fe , 



In the previous formula, the symbol = A:} denotes the number of components of the vector s 
equal to k. It follows that the term 2™=i' B a S a in Eq. ( |2.1| i is diagonal in the basis ( |2.1 l| i, with 
eigenvalues 

m-l m-1 m-1 

J] c a (s)B a = j] Z?„#{i,- = a} - #{*,■ = m) j] B„ . (2.13) 

a= 1 (2= 1 a=l 

In view of the latter equation, it is convenient to introduce the notation 

m-1 

B m = ~Y J B a , (2.14) 

so that Y!a=i Ba = 0. Using this notation we can rewrite the eigenvalue ( |2.13| l in the more 
compact form 

m N 

J] 5 a#{ii = ff} = 2 2? «- (2 ' 15) 

From the previous considerations and Eqs. |ZT2|-(|ZT5| it then follows that the Hamiltonian j2.l\ 
acts triangularly on the basis ( |2.1 1} , with eigenvalues 

2a I N 

E {tS = E e +--[j\t\-Y i B s \. (2.16) 

Using the previous equation, it is a straightforward matter to compute the partition function 
of the spin Calogero model ( |2.1[ ). To this end, let us represent the multiindex t in Eq. ( |2.17| > as 



(2.17) 




with A\ > A 2 > ■ ■ ■ > A r > and (k\ , . . . , k r ) = k e <P N . We then have 

Z(2aT/J) = q^Yj Z JklA ' Z «~ Z£ ' **' ' (2 ' 18) 

kePw ,li> ->i r JsO sef 

where the last sum is extended to all spin quantum numbers s compatible with the multiindex 
t given by Eq. ( |2.17) , i.e., such that s, ~>t whenever 1/ = Z, +I . In fact, since the latter sum 
clearly depends on ( only through k, from now on we denote it by S(k). Clearly, by Eq. ( |2.17| > 
we have 

£(k) = f] J] q r . (2.19) 

(=1 S] ^■■■■~^Sk j 

With this notation, Eq. ( |2.18| > becomes 

keP N Ai>->A r >0 



Z(2aT/J) = q^Yj S(k) Z q& 



The inner sum in the RHS is easily computed by performing the change of indices v, = A, - Ai+\, 
where i = 1, . . . , r and A r+ \ = 0, since the new indices independently range from 1 - 6, r to oo. 
We thus obtain 

r-l jk 

Z * 2 '"-"-'/"^ll,!-- < 2 - 2 °) 

A,>->A r ^0 i=l H 

and therefore 

T 

Z(2aT/J) -q ^ S(k) «i _ , (2.21) 

where the positive integers ^T,- were defined in Eq. < |1 . 10| . From the freezing trick relation ( |2.8| l 
and Eqs. ( |2.9) -( |2~2"T| >, we immediately obtain the following explicit expression for the partition 
function of the PF chain ( |l.l| i-( [L7b] i in the presence of an arbitrary su(m) magnetic field: 

N-r 

Z(q) = J] S(k) q^ ]K > Y]( 1 - q JK ') . (2.22) 
keP N i=l 

The previous argument can be applied with minor modifications to both the HS and FI chains, 
thereby obtaining the following general expression for the partition function of the HS-type 
chains (fTTT)-(fr7Ti: 

N-r 

Z(q) = Yj S ( k > V LrJ mK,) F[ (1 - ^''^ ■ (2 ' 23) 
keP N i=l 

Note that when the magnetic field vanishes, by Eq. ( |2.19) 2(k) becomes ]~[/=i \k~)' anc ^ tne 
previous expression for X reduces to its zero field version ( [1.9) . 

A few remarks on the function Z(k) are now in order. Let us first note that the latter function 
can be expressed as 



X(k) = f] (2.24) 



(=1 
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(2.25) 



where 

o-(k) = 2^ 1 

The function cr(fe) can be considerably simplified by noting that the summation indices s\, . . . , Sk 
can be expressed as 



m, ...,m,m— 1 , . . . , m — 1,.. 



(2.26) 



with 5£L, p a = jfc. Thus 



and therefore 



n = l 



The latter sum can be easily evaluated by noting that 

oik) = h k {q- B \.. 



(2.27) 



(2.28) 



where %(x) denotes the complete homogeneous symmetric polynomial of degree k in m variables 
...,x m ) = %, given by 

h k (x) = ^ *f---4 m . (2.29) 

Recalling that hk is the Schur polynomial associated with the single-row partition (k) and using 
Jacobi's definition of the latter polynomials in terms of determinants [42 1 we conclude that 





■ q- 


q -(m-2)B, 


■• q 


q -(m-3)B, 


q 






q~ B ' 




1 





7 -(fc+m-l)fl„ 
-(m-2)B,„ 

I 

7 -(m-3)B,„ 



(2.30) 



Expanding the determinant along the first row we obtain the alternative expression 

m m 

<T(k) = Y jq -^ B 'Y\( q - B '- q - B r 1 . 



(2.31) 



In particular, for spin 1/2 (m = 2) we have B\ — B — -B2, and either expression for cr(k) 
easily yields 



cr(k) 



q W)B _ q -(k + l)B ^ sinh ( (jt + 1)jgB ) 

<7 B - <7~ B 



(2.32) 



smh(j3B) 

The RHS of the latter formula can be conveniently expressed using the ^-number notation. Re- 
call, to this end, that given two real numbers x and w > the symmetric w-number [[x]] w is 
define by 

l[x]] w = —, r , 



W 2 — W 2 

9 



which reduces to the ordinary number x for w — I. We then have (in the spin 1/2 case) 

<r(t) = [[t + l]] 9 », 

so that the partition function of the spin 1/2 chains ( | 1 . 1 [ >-( fT/7] > can be concisely written as 

Z(q) = ^ fl[[*i + 1]] " 2S ■ /£W) ft 1 " ^ • (233) 

keP N i=l i=l 

Note, again, that in the absence of a magnetic field the previous expression obviously reduces to 
Eq. {L9]l with m = 2. 

3. The associated vertex models 

In this section we shall prove that the spectrum of the HS chains ( |l.l| i-( fL7) i coincides with 
that of a classical inhomogeneous vertex model, which differs from the one in Eq. ( |1.12| i by the 
addition of a magnetic field term. Our approach is based on the notion of generating function for 
the partition function, as used in Ref. OBI for the zero magnetic field case. 

To this end, following the latter reference we define the generating function T of the zero- 
field partition function ( |1.9| > as 

/- N-r 

T(x) = J] f[ h kl (x) ■ q"- ^ Y\(l - q J ^) . (3.1) 

kef N i=l i=l 

Since 

A t (l,...,l) = (" M ^ 1 )' 

by Eq. ( |1.9| > we have Xoiq) = T(\, ■ ■ ■ , 1). More generally, substituting the expression ( |2.28| > for 
cr(k) into Eqs. ( |2.22| >-( [2~24] > and using the definition of ^"(x) we readily obtain the identity 

Z(q) = T(q- B, ,...,q- B >"). (3.2) 

Similarly, the generating function T w for the partition function of the classical vertex model ( |1.12[ ) 
was defined in Ref. 1 35 1 by 

m 

T v (x)= J] i; ,|(n) -<» w /» w , (3.3) 

where the nonnegative integers w a (n) are given by 

w a (n) = #{n k = a] . 

Again, it is obvious that the partition function Z,^ of the model ( |1.12| >-( f]"T3) ) is the value of its 
generating function at the point (1, . . . , 1). One of the fundamental results in Ref. l35l is the fact 
that the generating functions p.l[ ) and ( |3.3[ > are identically equal, i.e., 

r(x) = T v (x) , Vx e R' n . (3.4) 
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Evaluating the previous identity at the point (1, . . . , 1) we obtain the equality of the zero-field 
partition functions Zo an d Z% > which is in fact the main result in Ref. (35]. On the other hand, 
if we evaluate the identity ( |3.4) > at the point (q~ Bl , . . . , q~ B,n ) and use Eq. ( |3.2| ) we immediately 
obtain 

m 

Z(q) = T V (q- B <,...,q- B >")= J] Z ™=' ^V ^ • (3-5) 

«i,,..,n.jv=0 

Taking into account that 

m N 

^ B a w a (n) = ^ B n , , 

a=\ !=1 

we can rewrite Eq. ( |3.5| l as 

m 

Z(q)= £ ^W-2£.*. S Z V ( 9 ), (3.6) 

«1,...,«A'=0 

where 2 V denotes the partition function of the inhomogeneous classical vertex model with en- 
ergies 

N 

E(n)=£ (n) n* e {l,...,m}, (3.7) 

i=l 

with £ (n) given by Eqs. (Tl2) -l |lTl3 



Therefore, as stated at the beginning of this section, the 
spectrum of the spin chains ( |1.1| )-( |1.7 1 is identical to that of the classical vertex models defined 
by Eq. ([Ll2)-<[377). 

In the spin 1 /2 case, the previous equation can be simplified by noting that 5{j, k) in Eq. ( |1.13| l 
can be expressed as 

6<J,k) = (j-W-k), j,k=l,2, (3.8) 

and similarly 

Bj = (3-2j)B, j=l,2. (3.9) 
Introducing the spin variables crj = 3 - 2n^, Eq. ( |3.7[ ) becomes 

. N-l N 

E(o-) = -2£(0(l-^ ! )(l+cr i+1 )-B^ C r i , cr h e {±1} . (3.10) 

!=1 1 = 1 

The latter equation can be alternatively written as 

. N-l N 

E{cr) = ~jYa £ «(<W + i - 1) - Yu fl ® ai ' (3 ' ! 1} 

i=i i=i 

where 

S(i)=B+^[8(i)-8(i-l)] (3.12) 

and we have set £(0) = S(N) = 0. The last two equations define a classical inhomogeneous 
Ising model in one dimension, where the coupling between the spins i and i + 1 is proportional 
to the dispersion relation £(/), and the external magnetic field (also inhomogeneous) is given by 
Eq.<[3TT2). 
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4. Thermodynamics of the spin 1/2 chains 



We shall next take advantage of the representation ( |3.10[ ) of the spectrum of the HS-type 
chains ( |l.l| i-( fL7] i with spin 1/2 in order to determine their equilibrium thermodynamics in a 
unified way. To this end, we must first normalize the Hamiltonian ( |1.6| l to ensure that the energy 
per spin is finite (and nonzero) in the thermodynamic limit. It is well known that when N » 1 
the mean energy of the HS-type chains is 0(N r ), with r — 2 for the PF chain and r = 3 for the 
HS and FI chains [19 3311341 . Hence from now on we shall take 

( K/N 2 , for the HS and FI chains 
J = \ (4.1) 
1 K/N , for the PF chain, 

where the constant K is independent of N. With this proviso, Eq. ( |3.10| ) can be rewritten in terms 
of the discrete variables x,- = i/N e (0, 1) as 

K N ~ l N 
E(cr) = - ^ Bi(l ~ o-d{\ + ir w ) -B^o-i, (4.2) 

where 



4 

!=1 !=1 



Si = < 



x, (l - xi) , for the HS chain 

Xi , for the PF chain (4.3) 

XiijN + Xi) , for the FI chain, 



and y N = (a - X)/N. We shall further assume that has a finite limit y ^ as N — > oo. With 
the above notation, the partition function of the chains ( |1.6| )-( fL7] i can be collectively expressed 
as 

£(q) — q^^' ■°" 2 ^ El ~ f +q " 2 ^ ... qd(o~N-\,o-N)EN-\-j (o-«-i+o"«)^-f (o-i+o-jv) (4 4) 

where 

d(<r,oO= f(l-£r)(l+cr'). 



Equation ( |4.4| ) can be more concisely written as 

Z(q)=tr(UT 1 ---T N - 1 ), (4.5) 
where U and T, are 2x2 matrices with elements 



= q-i^ , (r,W =/^ £ '-f (cr+ ^; < r,(7' = ±l, (4.6) 



or equivalently, 



-if;.)' *-(??'■ 



The transfer matrix T, has two distinct eigenvalues 



A u± = cosh(j8B) ± ^sinh 2 (J3B) + e- K P B ' , (4.8) 
and is therefore diagonalizable. A straightforward calculation shows that 



T t = PiDiPJ 1 , (4.9) 
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where 



and we have set 



We thus have 



D, 



V 
Ai- 



P, 



—a + r,- —a — r, 
1 1 



= sinh(ySB) , n = ^sinhOSB) 2 + e-^ £ < . 



where C, = P. l Pi+i is a symmetric matrix explicitly given by 



The key observation at this point is that 



( n + r M n - r M \ 

n - n+i n + r M 



Q=I + OCAT 1 ) 

as N — > oo. Indeed, note first of all that 

^ £j ^ £max i 

where 

r . 

, for the HS chain 

1 , for the PF chain 

yjv + 1 , for the FI chain, 
so that e, = O(l). On the other hand, since < x, < 1 the difference 



- Sj = < 



5y (1 - X; - Xi+i) , for the HS chain 
i , for the PF chain 



I N 



1 (Tat + *i + x i+i) > f°i" me FI chain 



is 0(N and hence 



— = 1 + OCAT 1 ), 



which establishes our claim. 

From Eq. ( |4.13[ ), it follows that when N ^ oo the leading term of Z.{q) is given 



where 



Note also that 



\x(p- N \UP l D 1 --D N - X )=\i 

N-l 

a ± = Y] k 



A+ 
A- 



i=i 



lim s\ = 0, lim en-\ = £oo = 

A 7 — >oo >oo 
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o, 
1 , 



for the HS chain 
for the PF chain 



y + 1 , for the FI chain, 



so that the matrices P\ and P^-i both have finite limits as N — > oo. Calling 

= ^sinh 2 (BB) + e-^ £ ~ 

we then have 



lim P^^I/Pi = 



1 

2r« 



1 r m + a 
-1 roo - a 



e"^ 1 
1 ^ 



1 1 



cosh(y3B) ( r M + cosh(£B) 
- coshOSB) 



(4.15) 



and therefore 



tr 



/l + 
A- 



= coshOSB) 1 + 



shQ3g) \ 



is the leading term of Z(q) as Af — > oo. Since the first two factors in the RHS of this formula are 
independent of N, in the thermodynamic limit the free energy per spin f(B, T) is simply given 
by 



1 7 1 N ~ 1 
f(B, T) = - lim = — lim AT 1 V logi ! + 

/= l 

1 N_i I 

= — lim AT 1 V log(cosh(0fl) + Jsinh 2 0SB) + . 

f3 n^oo — ^ \ y 



). 



(4.16) 



where the dispersion function s(x) is defined by 

x( 1 - x) , 
£(x) = < x , 

x(y + x) , 



for the HS chain 
for the PF chain 
for the FI chain 



(4.17) 



(cf. Eq. ( |4.3[ Q. We thus obtain the following remarkable formula for the free energy per site of 
the HS-type chains ([T76T>-(fTT7]>: 

f(B,T) = -- f log(coshOSB) + Jsinh 2 (BB) + ) dx. (4.18) 

P Jo 



From the latter equation it is immediate to obtain the magnetization per site: 



KB, T) = - 



df(B, T) 
dB 



■Jf 

Jo 



sinh(£B) 



^sinh 2 ^) 



: dx . 



(4.19) 



+ e -K/3e(x) 



For finite temperature, the magnetic susceptibility is also easily computed from the previous 
formula: 

dp.(B, T) f 1 cosh(BB)e-^ eW 



X(B,T) 



dB 



f 

Jo 



(sinh 2 0SB) + e-W*)) 



3/2 



dx . 



(4.20) 



In the case of the FI chain, some of the above equations have previously appeared in Ref. (6). 
More precisely, Eq. ( |4. 1 8] > with B = reduces to Eq. (21) of the latter reference, while Eq. ( |4.19| l 
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for K > essentially coincides with Eq. (24) of Ref. [6|, which is stated without proof. As 
first noted in this reference for the B = case, Eqs. ( |4. 18| l-( |4.20| i imply that the HS chain is 
thermodynamically equivalent to an FI chain with y = -1 and K — > -K, 

The free energy and magnetization per spin ( |4.18| l-( |4~T9| ) admit a simpler expression in terms 
of the so-called dressed energy Sfc(x) defined by 



namely, 



KB , . 

sinh(y6e dr (x)) = e^*'' sinhl/Jfi) , (4.21) 



i 



fi(B, T) = I tanh^Sd,-^)) dx . (4.23) 



These are the well-known expressions for the free energy (up to an additive constant depending 
on the normalization of *H) and magnetization per particle of a system of free bosons with one- 
particle energies e(x), x e [0, 1], in the presence of a constant magnetic field B. Note in this 
respect that, in contrast with the dispersion relation e(x), the dressed energy Sd r (x) depends on B, 
B and K. Note also that in the case of the HS chain Eq. ( |4.22[ ) is essentially equivalent to Eq. (37) 
in Ref. 11431 . deduced by means of Haldane's spinon gas approach ll23l . 

For future reference, we shall also compute the internal energy and entropy per spin. To 
begin with, the internal energy per spin u is given by 

-«n->wn)-i ff -^^L -M^A, ( 4,4) 

oB x 2 Jo \cosh(BB) + r(x) ) r(x) 

with 

r(x) = -yJsmh 2 (BB) + e^M . 
The entropy per spin s is then easily derived from the usual formula 

u(B, T) - f(B, T) 
s(B, T) = t JK . (4.25) 

4.1. The PF chain 

We have not been able to express the thermodynamic functions ( |4.18| >-( |4.25| l in terms of 
known special functions for arbitrary values of B and T, with the remarkable exception of the 
Polychronakos-Frahm chain. Indeed, we shall now show that in this case the free and the in- 
ternal energy, the entropy (and the specific heat) can be expressed in terms of the dilogarithm 
function EH 

Li 2 (z) = - f l0g(1 ~ ?) dt, zeC\(l,co), (4.26) 
Jo t 

while the magnetization and susceptibility are elementary functions of B and T. In Eq. ( |4.26[ ) log 
denotes the principal value of the logarithm (i.e., —n < Imlogz < n), and the integral is taken 
along any path joining the origin with the point z that does not intersect the branch cut (1, oo). 
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We start by performing the change of variable f 2 = a 2 + e **, where a - sinh(/JB) and k = KB, 
in Eq. ( |4.18[ ) for the free energy. We readily obtain 



f{B,T)= f (-L + J_)log(r + c)d*, 
J c \t + a t—aj 



>V5 

KB 2 



where c = cosh(J3B). The linear change of variable 

z = — = -e ±/3B (f + a) 

c + a 

makes it possible to evaluate each of the last two integrals in closed form: 

rV^ i og ( f + c ) ^ = p^ (±fl+ V^) / lQg(1 -z) t BB\ (J _ 
Jc t±a J- e ±2HB \ z z J 

= Li 2 (-e ±2BB ) - Li 2 ( - e ±/3B (±a + ^a 2 + e"*)) 
+ BBlog(±a + Va 2 + e"*) + B 2 B 2 . 
Taking into account the dilogarithm identity [44] 

2 2 

Li 2 (-e*) + UaC-e"*) = ~ - - , (4.27) 

after some straightforward algebra we finally obtain the following closed formula for the free 
energy of the PF chain: 



-KB 2 (f(B, T) + B) = — + 2ySBlog(sinh(j8B) + ^sinh 2 (BB) + er K P ) 
+ Li 2 (-e /JB [sinhOSB)+ ^sinh 2 (BB) + e~ K P ]) 
+ Li 2 (-e^ B [ - sinh(/3B) + ^Jsinh 2 (BB) + e"*?]). (4.28) 

The previous formula generalizes the analogous expression for zero magnetic field in Eq. (23) 
of Ref. [5|. Equation ( |4.28| l can be somewhat simplified by expressing the RHS in terms of the 
dressed energy. Indeed, by Eq. ( |4.21| > with s(x) — x we have 



ih(BB) + yjsinh z (BB) + er K P = ^(-f ±e *(D) , (4.29) 



and therefore 

-Kff{B, T)=j+ 2B 2 Bs dl (l) + Li^-e^^H)) + Li 2 (-e^( B+e * (1)+ f )) . 

From the usual formulas for the internal energy, the entropy and the specific heat (per site) 

du(B, T) 



c(B,T) 



dT 



in terms of the free energy, it immediately follows that these functions necessarily include a 
term proportional to the last two lines in Eq. ( |4.28| l, and therefore cannot be expressed purely 
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in terms of elementary functions. On the other hand, since the magnetization is obtained from 
the free energy by differentiating with respect to the magnetic field B, it is clear from Eq. ( |4.28| l 
that fi is an elementary function of B and T. Differentiating Eq. (438} with respect to B, or 
simply performing the change of variables t 2 - a 2 + eT K P x in Eq. ( |4.19| l, after a straightforward 
calculation we obtain the remarkable formula 



u(B, T) = 1 + — 
FV KB 



log(sinh06B) + ^sinh 2 (BB) + erW) - BB 



Using Eq. ( |4.29[ ), we can rewrite the previous formula in terms of the dressed energy as 

fi(B,T) = ^(s iT (l)-B). 



(4.30) 



(4.31) 



Differentiating the previous expressions for \i we also obtain the following closed formulas for 
the susceptibility per site: 



X(B,T) = 



K 



cosh(J3B) 



g 1 

■^jsmh 2 (J3B) + er K P 



2 / tanhQ3edr(l)) 
K \ tanh(ySB) 



- 1 



(4.32) 



5. The zero field case 

It is of interest to study the limit of the thermodynamic functions computed in the previous 
section for both B — > and T — > 0. We shall deal in this section with the first of these limits, 
leaving the analysis of the zero temperature limit for the next one. 

We shall start with the free energy and the internal energy per site, which can be immediately 
obtained by setting B = in Eqs. ( |4.18| l and ( |4.24| i: 



f(°> r > = -jjf io g (i+ e -? eW )d*, 
k r l s 

2 Jo 1 + , 



u(0, T) 



• e 2 



(5.1) 
(5.2) 



From the above formulas and Eq. ( |4.25| l it follows that the entropy per spin is given by 



s(0,T) 



KBs(x) 



Ax 



(5.3) 



= f log(l + e-^)) 1 

= JT 1 {log[2cosh(f s{x))] - ^ tanh(f eW )}dx. 

In the case of the HS chain, the latter equation is essentially equivalent to Eq. (15) in Ref. 
which was deduced using the equivalence of the HS chain to an ideal gas of spinons. The zero- 
field specific heat 

du(0, T) 

c(0, T) = 



dT 

can also be easily computed from Eq. ( |5.2| i, with the result 

^=(f) 2 ^ e 2 Wsech 2 (f eW )d,. 
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(5.4) 



Obviously, the zero-field magnetization vanishes (cf. ( |4.19[ )), while Eq. ( |4.20| i for the suscepti- 
bility reduces to 



X(0,T)=/3 



e- e(x) dx. 



(5.5) 



In fact, the zero-field susceptibility can be expressed in terms of standard special functions for all 
spin chains of HS type. For the PF chain, we have already obtained a formula for^- in terms of 
elementary functions for all values of B and T, cf. Eq. ( |4.32| i. For the HS chain, a straightforward 
calculation shows that 



X(0,T) = 



2rB 



e *W8 x , 



erf(^f), K^O 
erfi^), K i0 , 



(5.6) 



where 



erf (z) = 



dt , erfi(z) = 



e dt = -ierf(iz) 



denote respectively the error and the imaginary error functions ll45l . Similarly, for the FI chain 
we have 



2\K\ 



vfi((y + 2)Jf )-erfi( rA /f ), K>0 
erf( (r + 2 )A /T)-erf( rA /^), *<0. 



(5.7) 



From Eqs. ( |4.32| i, ( |5.6| ), and §5.7) , and the well-known asymptotic formula [45 
erfiz 



e~ 



it follows that when K > and T — > we have 



X(0,T) 



/ 2"jg p/f/3/8 



2 Kfi/2 
K C 



(y+2)K 



^■y+DKfl/2 



, |argz| < - , 

for the HS chain 
for the PF chain 
for the FI chain . 



Again, in the case of the HS chain the previous asymptotic formula was first obtained by Haldane 

in Ref. [23 ]. From the latter equation it is apparent that the ferromagnetic zero-field susceptibility 
(and hence the correlation length) is exponentially divergent as T — > 0. This shows that all 
ferromagnetic chains of HS type undergo a phase transition at T — B — 0, which we shall 
analyze in detail in the next section. 

Equations \53\ and \5A\ imply that the zero-field entropy and specific heat are even func- 
tions of the coupling constant K. This fact, which is not obvious a priori, was first noted by 
Haldane ll23l for the entropy of the original HS chain. Haldane also pointed out that (with an 
appropriate choice of the ground state energy, which amounts to replacing 1 - Sy by Sy in 
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Eq. ( |1.6| l) the internal energy of the HS chain is again even in K. One can easily show that this 
property is shared by all spin chains of HS type, since from Eq. \5.2\ it readily follows that 



with 



K Ke 
u(0, T; K) - h(0, T; -K) = — e(x) dx = — - , 

2 Jo 2 



so 



( . for the HS chain 

2 , for the PF chain (5.8) 



i (3y + 2) , for the FI chain. 



In other words, u(0, T) - is an even function of K, which can be seen to be equivalent to 
Haldane's statement given our normalization of Iri. A similar property, proved in exactly the 
same way, also holds for the free energy /(0, T). 

As we saw in the previous section, in the particular case of the Polychronakos-Frahm chain 
the thermodynamic functions ( |5.1[ )-( |5~4"| > can be expressed in closed form in terms of the diloga- 
rithm function. Indeed, setting B — in Eq. ( |4.28| l we obtain the simple formula 

/(0,r) = -^(Li 2 (-e-^ 2 )+^). (5.9) 

From the latter equation it is straightforward to derive similar expressions for the internal energy, 
the entropy and the specific heat of the PF chain in the absence of a magnetic field, namely 

«(0,D = JL(li 2 (-e-W 2 )+ g)_ I io g (i +e"^ 2 ) (5.10) 
V( °- n ±( L i 2 (-e-^ 2 )+ ^)-log(l + e-^ 2 ) (5.11) 



5.1. Discussion 

Using the previous explicit formulas, we have plotted the thermodynamic functions of the PF 
chain for zero magnetic field ( |5 .9[ >— ( |5TT2] > in terms of the dimensionless temperature r = (\K\By 1 ; 
see Fig. [T] 

It is immediately apparent the qualitative similarity of these plots with the corresponding 
ones for a two-level system like, e.g., a paramagnetic spin 1/2 ion or the one-dimensional Ising 
model with zero magnetic fielcj^][46|. In particular, the specific heat c(0, T) clearly exhibits the 
so-called Schottky anomaly, characteristic of this type of systems. The Schottky peak, which 
can be numerically computed without difficulty from Eq. ( |5.12) , is located at \K\B 0.127762, 
with a maximum value of c(0, T) approximately equal to 0.439229 &b- On the other hand, the 
low temperature behavior of s(0, T) and c(0, T) markedly differs from that of a two-level system. 



4 For the latter model, the effective half-gap is equal to the coupling between the spins. 
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Figure 1 : Free and internal energy, entropy, and specific heat (all of them per site) of the PF chain with zero magnetic field 
(solid blue line) and of a two-level system with half-gap A = \K\/S (dashed red line) as functions of the dimensionless 
temperature r = (\K\f3y 1 . In the first two plots, uq = - sgnA") denotes the common value of u(0, 0) and /(0, 0), 

cf. Eq. (5.17). Insets: low-temperature behavior of each of the above thermodynamic functions. 



Indeed, from Eqs. ( |5.1 l| l and ( |5.12| i it immediately follows that 

s(0, T) n 1 



c(0, T) _ n 2 
k B ~ 3\K\P 



+ 0(r z ), T->0, 
+ 0(T 2 ), T^Q. 



(5.13) 
(5.14) 



Hence both functions have a nonvanishing first derivative at T — 0. As is well known, the entropy 
s&{T) and the specific heat ca(T) of a two-level system with energy gap 2A > satisfy 



sa(T) 



= log(2 cosh(/?A)) -jSAtanhOSA) ~ 2/3 A e 



-2/3A 



T -> 0, 



(5.15) 



ca(T) 



08A) 2 sech 2 (j6A) ~ 4/3 2 A 2 e" 2/3A , T -> , (5.16) 



so that the derivatives of all orders of s&(T) and ca(T) vanish at T = 0. 

We shall next show that the previous analysis essentially holds for the HS and FI chains as 
well. To begin with, note that from Eqs. ( |5.1| l-( |5T2"| ) it follows that 



/(0, 0) = u(Q, 0) = (1 - sgn K) = u , 
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(5.17) 



and thus 



f(0,T)-uo = -- log l+e-- eW )dx, (5.18) 
P Jo 

.['- 

Jo 1 + ( 



m(0,7)-« = — sr— dx. (5.19) 

2 Jn 1 xk^sW 



On the other hand, the free energy f A and the internal energy u& of a two-level system with gap 
2 A, when normalized so that «a(0) = /a(0) = 0, are given by 

MT) = -i log(l + e^ A ) , ua(T) = — ^- . (5.20) 

From the latter equation and Eqs. ( |5.15| l-( |5.16[ ), it follows that the thermodynamic functions of 
a two-level system with gap 2A are obtained by taking s(x) in Eqs. ( |5.3| l-( |5~4] i and ( |5.18| l-( |5.19| l 
as the constant 

4A 

*- m . (5.21) 

This observation suggests that, if the ratio A/|*| is chosen appropriately, the thermodynamic 
functions of a two-level system with gap 2A and those of a spin chain of HS type with coupling 
* (and zero magnetic field) should behave in a qualitatively similar fashion. In order to establish 
the relation between A and *, we note that in the high-temperature limit one has 



lim [m(0, T) - Mnl = , lim ua(T) = A . 

T-tco L J 4 T-tco 

For both limits to coincide, we must have 



A_ _ so 
1*1" 4 ' 

i.e, £a should be equal to the average of s(x) over the interval [0,1]. In particular, in the case of 
the PF chain we have eo = 1/2 (cf. Eq. ( |5.8[ )), which explains the choice A = |AT|/8 in Fig.[T] 
As an example, we have plotted in Fig. [2] (left) the specific heat per site of the PF, HS, and FI 
chains (with y = 0, 1 in the latter case) compared to its counterpart for a two-level system with 
half-gap A = 1*1/8. Here * denotes the coupling of the PF chain and, following the previous 
analysis, we have fixed the corresponding couplings *, of the other chains by the requirement 
* J^ 1 Si(x) dx = * J Q l xdx = K/2. 

We shall next derive the low temperature behavior of the thermodynamic functions of the HS 
and FI chains in the absence of a magnetic field. Our strategy shall be to obtain an asymptotic 
formula for the specific heat ( |5.4[ ), which we shall then integrate to obtain similar formulas for 
the remaining thermodynamic functions. 

Consider first the HS chain. Since in this case s(x) — x(l - x) is symmetric about x = 1/2, 
we can write Eq. ( |5.4[ ) as 

C (°> T ) „,2 f' /2 s 2 (x)e-' lE(x) 



2,\ z f — — — — -dx. 
Jo 



(1 +& -Mx))2 
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where A = \K]/3/2 is a large parameter. Performing the change of variables t - A s(x) we imme- 
diately obtain 

c(0,T) 2 C m t 2 e-> dt 2 f°° t 2 e"' 

■ 0{A ) 



_ 2 r A /4 f 2 e~' df _2 r f 2 e" 
~IJo (1+e-') 2 J^Tf ~ ^ Jo (1+e-') 2 f + ' 



(see Appendix A for the details). The last integral is easily computed: 

Jo (1+e-') 2 ^ Jo ^ " 6 

(cf. Eq. ( |4.27| i), which combined with the previous equation yields 

c(0,T) _ 2n 2 2 



+ 0(T A ) . (5.22) 



Note that the leading term in this formula is twice that in Eq. ( |5.14[ ) for the PF chain, due to 
the fact that in the present case the leading contribution to the integral \5.A\ comes from both 
endpoints. Apart from this inessential difference, the specific heat at low temperature behaves as 
its counterpart for the PF chain, i.e., it increases linearly with the temperature. 
Consider next the FI chain. If y + 0, proceeding as before we obtain 

c(0,T) 2 f 1 e 2 (x)e-' lEW J _ 1 f 2 e~' df 



= A 1 ' 7 : ; — Ax = 



o (1 +c- '•''")- " A Jo (1+e-') 2 



1 P°° f 2 e-' n 2 
yAJ (1+e-') 2 3y\K\p 



Again, in the low temperature range the specific heat increases linearly with the temperature. 
Finally, for the FI chain with y = we have 

c(0,T) , f 1 x 4 e- vU? 1 C A f 3/2 e-' 

'- 1 df 



, r x'q-^ l r 

Jo (l+e-^ 2 ) 2 2 VI Jo 



1 f M f J / 2 e-' , 
VI Jo (1+e-') 7 



2 VI Jo (1+e-') 2 

f 3 / 2 e-' 
2 VI Jo d+e-') 2 

where the last integral can again be exactly evaluated: 

,3/2 -r 00 /-oo 

f 6 9 dr^V(-ir'« ^V'dr 
o (1+e-') 2 £f Jo 

- e X" f3/2 e ' df = ^ r ™ = ^ ^ • 



where 



(-1)" 



11=1 
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is Dirichlet's eta function. Using the elementary relation 



7 7 (z) = (l-2 1 - z )^(z), 
where f (z) is Riemann's zeta function, we finally obtain 



r 

Jo 



< 5 '-' 1 K(3/2) , 



(1+e-') 2 4 \ V2 

and hence 

= j(V2- l)?(3/2) S^ + 0(Pz- mi2 ). (5.24) 

In contrast with the previous cases, the specific heat c(0, T) has now an infinite first derivative 
at T — 0; see Fig. [2] (right). Hence the FI chain with y — B — exhibits a second-order phase 
transition at zero temperature. 

Since asymptotic expansions can in general be integrated termwise to yield valid asymptotic 
expansions [47 1, integrating the asymptotic formulas ( |5.14| i and ( |5.22| i-( |5.24| i we can readily 
derive the low temperature behavior of the remaining thermodynamic functions of all spin chains 
of HS type. Thus, for the FI chain with y — we have 

/(o, T) = uo - -j /r 3/2 + oG(rV W2 ) (5.25) 

m(0, T) = uo + y ^ + °0 3 " le "' W2 ) ( 5 - 26 ) 
= 2t,/T 1/2 + 0(e-'^ 2 ) , (5.27) 



where we have set 



tj = I ( V2 - 1) f (3/2) Jj^ - 0.719227 



1^1 



-1/2 



For all other chains of HS type, the analogous results can be concisely summarized as follows: 

2 

/(0,r) = M0 -^2-+O(r 3 ) (5.28) 



where now 



«(0,r) = iio + 7 j^+o(r J ) (5.29) 

+ 0(J 2 ) , (5.30) 



6|^ 2 

5(0, T) _ 777T 2 - 



2 , for the HS chain 
1 , for the PF chain 
l/y , for the FI chsiawitky + 0. 
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Figure 2: Left: specific heat per site of the PF (blue), HS (light blue), and FI (y = 1, orange; y = 0, green) chains with 
zero magnetic field and of a two-level system with half-gap A = \K\/S (dashed red line) as functions of the dimensionless 
temperature r = (|jf|/?) _1 . Right: specific heat per site of the FI chain with y = (solid blue line) versus its low- 
temperature approximation (5.24) (dashed red line). 



6. The zero temperature limit 

In this section we shall analyze the low temperature behavior of the thermodynamic functions 
of a spin chain of HS type in the presence of a nonzero magnetic field. In the ferromagnetic case 
(K > 0), we already know from the previous section that there must be a phase transition at 
T — B — 0. This fact is also apparent from Eq. ( |4.19[ ) for the magnetization, which for K > 
implies that 

limyu(fl,r) = sgnfi. (6.1) 

The latter result can be corroborated by a computation of the zero-temperature magnetization per 
site directly from its definition. Indeed, from Eqs. (|4.2|i-(|4~3"|l it follows that 



KB,0)= lim i/jVA , (6.2) 
N— >oo Jy \r- ' / 

,= 1 g 

where ( ■ ) g denotes the average over all ground states. When B + and K > the ground state 
is obtained by taking cr,- = sgn B for all i in Eq. ( |3.10| >, and is therefore nondegenerate. Hence 
for B + we have 

fi(B,0)= lim — • Nsgnfi = sgnfl, 

N->co 

as before. On the other hand, when B — the ground states are obtained from sequences of the 
form cr = (1, , . . , 1, — 1, . . . , — 1). The average of Yjf = i 0"j over such sequences clearly vanishes, 
since the contributions of the sequences with k 1 's and N—kVs cancel each other. Thus /i(0, 0) = 
0, again in agreement with Eq. \6.\\ . 

In the antiferromagnetic case, it can be shown from Eq. ( |4. 19| l that there is a phase transition 
at zero temperature and magnetic field B = +B S , where the saturation field B s is given by 



B s = i \K\ max s(x) = 

2 *e[0,l] 



for the HS chain 



If , for the PF chain (6.3) 

If (y + 1) , for the FI chain. 
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Indeed, when K < we can rewrite Eq. ( |4. 19| l as 

MB, n = kJ 1 " 1 [i + 4(1 - e -w r v^]~ 1/2 dx t 

where A = \K\B — > oo, 

1 1 , for the PF and FI chains 

K = \ 

2 , for the HS chain, 



(6.4) 



and 

2\B\ 

<p(x) = s(x) - — = s(x) - b . (6.5) 
1*1 

Note that the value k — 2 for the HS chain is due to the symmetry of e(x) = x(\ - x) about 
x = 1/2 in this case. If \B\ > B s then b > no 

^j(x) < and e^ w -> as /I -> oo. Hence 



x = 1/2 in this case. If |Z?| > Z? s then Z? > max s(x) = max e(x), so that for all x e [0, k] we have 

xe[0,l] .te[0,^] 



p\lK 

\B\ > B, => lim |u(.B, r)| = /c dx = 1 , (6.6) 

Jo 

i.e, for \B\ > Z? s the magnetization per site saturates. On the other hand, if \B\ ^ B s we have 
(p(x) < for ^ x < xq and (p(x) > for xq < x ^ k, where xq denotes the unique root in [0, a:] 
of the equation tp(x) = 0, i.e, 

s(x Q ) = b , with x Q e [0, k] . (6.7) 

Thus in this case we have 

|0 , for < x < x Q 
oo , for Xq < x k , 



lim e^ (A) = 

A— >oo 



and therefore 

\B\ < £ s => lim MB, T)\=k dx = kxq . (6.8) 



Note that the RHS in both Eqs. ( |6.6[ ) and ( |6.8| > tends to 1 as \B\ — > B s , since clearly xq tends to 
1/*: as |B| -> fi s . 

For the HS chain we have e(x) - x{\ - x), and hence 



1 1 1 1 2\B\ 

X « = ~2-i~A~ b =~2-fA-W\- ^ 



Equations ( |6.3| > and ( |6.6| >-( |6T8) l (with k — 2) immediately yield 

sgns(l- w)' 



lim u(B, T) 

T^0 



(6.10) 

sgnfi, |fl|>f, 
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which was first derived by Haldane using the spinon gas formalism [23]. Likewise, for the PF 
chain we have k = 1, s(x) = x and thus xq = b = 2|B|/|^T|, so that 



'2B ini < ffl 



limju(fi,r) = f (6.11) 
\sgnfi, ||||>Jf. 

Finally, in the case of the FI chain k = 1, e(x) = x(y + x), 



and thus 



lim u(B, n = 



sgna^+f -y), |B|<lfl(y + l) 



sgnfi, |B|>f(y+l), 



in agreement with the result in Ref . (6) . 

From Eqs. ( |6.10| l-( |6.13| l it is apparent that, although the magnetization per site is continuous 
at B — +B S , the susceptibility^ = |g has a discontinuity at these points (cf. Fig.[5jl. Hence, in 
the antiferromagnetic case all three chains of HS type present a second-order phase transition at 
T — and B = +B S , where the saturation field Z? s is explicitly given in Eq. \63\ . The precise 
behavior of the zero-temperature magnetization and susceptibility for \B\ ^ B s is, however, quite 
different for each type of chain. Indeed, in the case of the HS chain the susceptibility is an 
increasing function of B for ^ B < B s , and diverges as (B s - |B|) -1 ' 2 when B — > ±B S +. On the 
other hand, for both the PF and FI chains the susceptibility has only a jump discontinuity at the 
critical points B = +B S . When ^ B < B s the susceptibility is constant for the PF chain, while 
for the FI chain it decreases monotonically. Finally, it should be noted that the first derivative of 
the zero-temperature susceptibility has a jump discontinuity at B — for the FI chain with y > 
and the HS chain, while it diverges as |B| -1 ^ 2 at this point for the FI chain with y — (cf. Fig . [3J . 

As in the ferromagnetic case, the zero-temperature magnetization can also be computed di- 
rectly taking advantage of the exact knowledge of the spectrum. Indeed, consider first the PF 
and FI chains. In this case the term proportional to K = -\K\ in Eqs. ( |4.2| i-( |4~3"j ) is minimized 
by the sequence <tq - (. . ., — 1, 1, . . . , — 1, 1, — 1, 1), characteristic of an antiferromagnetic sys- 
tem, while the term proportional to B is again minimized by the sequence with all components 
equal to sgnB. It is therefore clear that the ground state is obtained by flipping the compo- 
nents of (To equal to - sgnB starting from the left, until the resulting decrease 2\B\ in the term 
proportional to B in Eq. (|4.2[) is offset by the increase \K\s, = \K\e{xi) in the other term. Thus 



if \B\ > j\K\ max e(x,) = B s the ground state is again obtained by taking <r ; = sgnB for all 

i = l,...,N, and the magnetization attains its saturation value sgnB. On the other hand, when 
|B] ^ B s the ground state is obtained from the sequence 



(sgnB sgnS,-l, 1,-1,1,..., -1,1), (6.14) 

where for large the index i is (approximately) determined by the condition 



\K\sixd - 2\B\ 
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(6.15) 




By Eq. ( |6.2| i, the magnetization per site is therefore giverjj 

i 

ju(B,0) = sgnZ? • lim — = sgnB ■ lim x, . (6.16) 

N—>oo ]\f N—>co 

Taking into account that lim x, = xq, where as before Xq is determined by Eq. (16. 7b (with k = 1), 
we obtain the following formula for the zero-temperature magnetization per site for \B\ ^ B s : 

jt/(B, 0) = sgn B ■ xq . 

The latter equation clearly agrees with Eq. ( |6.8| l with k — 1 . For the HS chain the analysis is very 
similar, except that the sequence <Tq minimizing the term proportional to K in Eq. ( |4.2| > is now of 
the form 

0-q = (...,-1,1,..., -1,1,-1, 1,-1,1 ...,-1,1,...). 

[N/2] 

This is due to the fact that in this case the dispersion relation e(x) has a maximum at x = 1/2, is 
symmetric about this point, and is monotonically increasing in the interval [0, 1 /2]. For the same 



More rigorously, for large N there could be a small number of ground states with corresponding sequences almost 
identical to j6.14} , so that the average of 0"i over these ground states is approximately equal to its value on the latter 
state. 
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reason, if \B\ < B s the ground state is now obtained from a sequence of thefornf] 



(sgnfi, . . . , sgnB, —1,1, 



■1,1, sgnfi, . . . , sgnB) , 



(6.17) 



where i is again approximately determined by Eq. ( |6.15| l. Thus in this case Eq. ( |6.16| l should be 
replaced by 

2i 

yu(B, 0) = sgn B lim — - 2 sgn B lim x, = 2 sgn B ■ xo , 

which is essentially Eq. ( |6.8| l for k — 2. 

We shall next analyze the low-temperature behavior of the thermodynamic functions of the 
antiferromagnetic chains in the physically more interesting regime < \B\ < B s . To this end, we 
start by rewriting Eq. ( |4. 18| > for the free energy per site as 



-B(f(B,T) + \B\)= f log 
Jo 



i +e -wn 1 



dx, 



(6.18) 



where A = \K\B — > oo and ip{x) is given by Eq. ( |6.5[ ). It is straightforward to show that the 
function 



g(x, z) = log 



1±£ + I V(l -z) 2 +4e*W 



(6.19) 



satisfies 



az 2 



Integrating this inequality with respect to z over the interval [0, e 2iS ' s '] we obtain 

< g(x, e^ 1 ) - g(x, 0) < \ e-i^W e -^ |B| = I e^V^ 1 < ~ e^ 1 . 

Integrating now with respect to x over [0, 1] and taking into account Eqs. ( |6. 1 8| )-( |6Tl~9] > we easily 
arrive at the asymptotic relation 



-J3(f(B,T) + \B\) = log 



-I 



1 (l + Vl +4e^W A 



!/* f i + VI +4e*> ^ 



* log 



dx + 0(e^ |B| ) 
dx + 0(e~ m ) , 



(6.20) 



where we have used the notation introduced in Eq. ( |6.4) . Since i^(x) is negative for ^ x < xq 
and positive for xq < x ^ k, it is convenient to rewrite the latter equation as 



^ (f(B, T) + \B\ + K ^-£ <p(x) dx) = JT ° log 



f 1 + Vl +4e*) 



dx 



log 

J An 



2 V 4 



dx + 0(e~ m ) 



(6.21) 



As before, in general there could be a small number of ground states corresponding to sequences nearly equal to 
Eq. (6.17) , but this does not affect the conclusion for the reason given in the previous footnote. 
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Consider now the first integral in the RHS of the previous equation. Performing the change 
of variable t = -A (fix) and proceeding as in the previous section we readily obtain 



f 

Jo 



log 



( 1 + Vl +4e^W 



dx = 



1 



.,1/; 



( 1 + Vl + 4e- 



i r 01 

/le'(xo) Jo 



log 



log 

( 1 + Vl +4e-' 1 



df 



df + 0(A -2 ) , 



(6.22) 



where the error term has been estimated as in Appendix A Note that if < |Z?| < B s we must 
have < xq < k, and consequently s'(xo) > 0. Observe also that, in order for the above estimate 
to be valid, we must have Ab » 1, or (disregarding an inessential factor of 2) 

B\B\ » 1 . 

The second integral in Eq. ( |6.21[ ) can be similarly dealt with by the change of variable t - Aip(x): 



J 1 



log 



2 6 2 



<p(x) 



f Jl + -e-'W 



dx 



41 * 

<te (-*o) Jo 



df 



s'(x(t)) 



df + 0(/T 2 ) . 



(6.23) 



This formula is correct provided that /ty>(l/*c) = 8QK\s(1/k) - 2\B\) » 1, or equivalently (since 
\K\eil/K) = 2B B ), 

B(B S - \B\) » 1 . 

Combining Eq. ( |6.21| i with Eqs. ( |6.22| i-( |6~2"3] > we obtain 



where 



-i" 

Jo 



log 



- (/(S, 7") + |5| + — f cp(x) dx) 
C 1 + Vl + 4e" ; 



+ 0(A- 2 ) , 



df, h 



-f 

Jo 



log 



-e- ,/2 + Jl+fe 



df 



(6.24) 



(6.25) 



and we have omitted an exponentially small term 0(e <§c 0(/l ). Using the value of these 
integrals in Eqs. ( |B.2| >-( |B~4] > of Appendix B we finally obtain the asymptotic formula 

„~2 



f(B, T) = 



k \k\ r l,K 

"2-1 ' 



ip(x) dx ■ 



6\K\B 2 s'(x ) 



+ 0(r J ) , 



provided that 



(6.26) 



(6.27) 



Since Eq. ( |6.26[ ) is an asymptotic power series expansion in B, it may be differentiated 
termwise P7l to yield corresponding expansions for the remaining thermodynamic quantities 

29 



as T — > 0. In the first place, taking into account that Xq and ip(x) are independent of the tempera- 
ture we immediately obtain 



K\K\ r l/ " KIT 2 

u(B,T) = -\B\-4-l ip(x)dx+ +0(T 3 ), (6.28) 

2 J X( , 6\K\B l e'{x {) ) 



'-to 

c(B, T) 



+ 0{T Z ) , (6.29) 



s(B,T) K7T 2 n 

k B 3\K\Bs'(xo) 

In particular, in the case of the HS chain Eq. ( |6.30[ ) is in agreement with Haldane's original 
result 1231 

In order to compute the asymptotic expansion of the magnetization, we observe that, by 
Eq. (|6T7), 

dxo 2 sgn B 
~dB ~ \K\s'(x ) ' 

and, since <p(xo) = 0, 

d r l/K _ C 1/K d(p(x) , 2 sgn B / 1 



We thus have 



Note, in particular, that the leading term in the latter expansion coincides with Eq. ( |6.8| >. Using 
the values of Xq and k in Eqs. ( |6.9| l-( |6~T2l i and ( |6.4| ) we easily obtain 



MB, T)\ =kx q - K fl '^ \ + 0(T 3 ) . (6.31) 



MB, T)\ = l- + (x _ M)- 3/2 + 0( r 3 ) , (6.32) 

for the HS chain, and 

1^, r)| = \ - | - ^ (r 2 + ftf 3 " + o(7*) , (6.33) 

for the FI chain. In the case of the PF chain the 0(r 2 ) term in Eq. ( |6.31| l vanishes. In fact, in this 
case Eq. ( |4.30| l more directly yields 

= 2JB\ _ cr^ + Jl e ^ M \ (634) 
F \K\ \K\B \B ) 

The approximations ( |6.32[ )-( |6.34| i are in excellent agreement with the exact result ( |4.19| l (or ( |4.30| l, 
for the PF chain) provided that the magnetic field is in the range ( |6.27| i; see, e.g., Fig. |4]for the 
HS and FI chains. 
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Figure 4: Left: absolute value of the relative error of the approximation (6.32) for HS the chain for different values of 
the dimensionless temperature. Right: similar plot for the approximation (6.33) for the FI chain with y = 1. 

7. Connection with the one-dimensional Ising model 

We have seen in Section [3] that the chains ( |1.6| l-( [L7] i with spin 1 /2 are isospectral to the 
classical inhomogeneous Ising model in an external (also inhomogeneous) magnetic field defined 
by Eqs. ( |3.11| i and ( |3. 12) >. In this section we shall exploit this fact to study the connection 
between the above chains and the standard (homogeneous) Ising model, whose Hamiltonian we 
shall conveniently take as 

N N 

'Hj = jYf\-SiS M )-BY l Si, S N+l =S x . (7.1) 

With this normalization, the free energy per site in the thermodynamic limit is given by [48 1 

fj(T, B) = - i log(cosh(£B) + ^sinhOSB)^ + ) , (7.2) 

which is clearly reminiscent of the analogous Eq. ( |4.18[ ) for the spin chains of Haldane-Shastry 
type. In fact, we can give a heuristic derivation of the latter equation from Eqs. ( |3 . 1 l| l-( |3~T2"] i and 
the previous formula for the free energy per site of the Ising model. To this end, note first of all 
that the latter equations can be written as 

i=l i=l 

where 

We saw in Section|4]that when » 1 the term e, - e,_i is 0(N~ l ), so that in this limit we have 

K N-1 N 

i=i i=i 
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Proceeding as in the latter section, from the previous equation we can show that when N » 1 the 
partition function of the chains ( |1.6[ )-( fL7] i is approximately given by 



JV-l 

where A+(J) denotes the largest eigenvalue of the transfer matrix the Ising model ( |7.1| > with 
coupling J. Hence, in the thermodynamic limit the free energy per site of the chains ( |1.6| >-( fL7| > 
is given by 

1 7 1 

f(B, T)=- lim ^± = - lim — V log A + (^) . (7.3) 

A<->oo _/Vp A<->oo _/Vp 

On the other hand, when N — > oo the partition function Xj of the Ising model is related to the 
largest eigenvalue A+(J) of its transfer matrix by |48 1 

Zj - A + (J) N , 

so that 

MB, T) = - lim 1^ = -i log ^ + (7) . 
Inserting this identity (with / = ^p) into Eq. ( |7.3| l we obtain 

i r 1 

f{B, T) = lim - V / & (2?, r) = /aa (B, r) d* , (7.4) 

JV->00 A/ Z— i 4 Jq 4 

where we have used the fact that e, = e(x,) = e(i/N). Equation ( |4.18| l for the free energy of the 
spin chains ( |1.6| l-( [L7| ) readily follows from the previous equation and the formula f7.2| i for //. 

Equation ( |7.4[ ) for the free energy of the Haldane-Shastry spin chains admits a statistical 
interpretation that we shall discuss next. Indeed, let us first rewrite the latter equation as 

f(B,T)= - f f &A (B,T)dx, 



K JO 



where k is defined in Eq. ( |6.4[ ). Since the function e(x) is monotonically increasing in the interval 

4 



[0,k], performing the change of variables J = in the latter integral we immediately obtain 



r»max(0,7 ) 

f(B,T)= p(J)MB,T)dJ, (7.5) 

«7min(0,7o) 



where 



K 



f| , for the HS chain 

f , for the PF chain (7.6) 



and 



f(y+l), for the FI chain , 
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Note that Jq has the sign of K, so that the integration range in Eq. ( |7.5| l is [0, Jq] in the ferromag- 
netic case and [7o, 0] in the antiferromagnetic one. Furthermore, from Eq. ( |7.7| i it immediately 
follows that 

-*max(0,/o) 



p(J)dJ= - I d/t= 1 

Jmin(0,y ) * Jo 



/min(0,7o) 

Hence, by Eq. ( |7.5| >, the free energy of each of the spin chains ( |1.6| l-( fL7| > is the average of the 
free energy of an ensemble of standard Ising models weighted by the function p(7) in Eq. ( |7.7| i. 
This function, which depends on the chain under consideration through the dispersion relation 
e(x), is easily computed. Indeed, for the HS chain we have 

1 / 1\ 2 47 
e(x) = x(l - x) = - - \x- -j = — , 

so that 

s'(x) = 1 — 2x — \ 1 . 

V K 

From Eq. ( |7.7[ > with k — 1/2 we easily obtain 

p(J) = . 2 = ; 1 ■ (7.8) 

JK(^-J) 2V7o(7o-7) 

Similarly, in the case of the PF chain s'(x) = 1 and therefore 



1*1 \Jo 

is constant. Finally, in the case of the FI chain we have 



4 1 

p(J) = = m (7.9) 



2 „,2 



e(x) = x(y + x) = (x+ |j - 



T 2 _ 47 
4~ ~ Y 



and therefore 



, , 167 
e (x) = 2x + y = + — 



Hence the weight function p( J) is given in this case by 

4 y + 1 

y/K(y 2 K + 167) ^ + 4(r + 1)-/) 

Introducing the dimensionless variables |7o|p(7) and j = J /Jo (where j € [0, 1] in the ferromag- 
netic case and j e [—1,0] in the antiferromagnetic one) we can rewrite the previous formulas 

as 

±(1-/T 1/2 , for the HS chain 

\Jo\p(J) = \l, for the PF chain (7.11) 

(T + l)[r 2 + 4(r + 1)/T 1/2 , for the FI chain . 
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Figure 5: Normalized weight function \ Jo\p(J) versus dimensionless Ising coupling j = J/Jq for the HS (light blue), PF 
(blue) and FI (y = 0, green; y = 1, orange) chains in the ferromagnetic case. 

Equation \1.5\ clearly implies that any thermodynamic quantity of a spin chain of HS type is 
the average of this quantity over an ensemble of standard Ising models with respect to the appro- 
priate weight function p( J) in Eqs. ( |7.8| l-( |7TT0] i. In other words, a single spin chain of HS type is 
thermodynamically equivalent to a suitably weighted ensemble of standard Ising models. From 
this novel point of view, the fundamental difference between the three types of HS chains \\.6\ - 
( |1.7| i is the fact that the weight function p( J) increases with the Ising coupling constant J for the 
HS chain, is constant for the PF chain, and decreases with J for the FI chain; cf. Fig. [5] The 
study of this weight function can also uncover some unexpected relations between the different 
families of chains of HS type. For instance, from Eq. \1 .1 \\ it follows that the normalized densi- 
ties of the HS chain and the FI chain with y — are dual, in the sense that they are related by the 
reflection / h 1 — j (cf. Fig.|5]l. 

8. Conclusions and outlook 

In this paper we derive the thermodynamics of the three families of spin chains of Haldane- 
Shastry type ( |1.6| l-( fL7] i in a unified way. The main idea behind our approach is to exploit the 
equivalence of these chains with a corresponding family of Ising-like inhomogeneous vertex 
models. In the absence of an external magnetic field, this fact had already been conjectured by 
Frahm |5| and by Frahm and Inozemtsev [6|, and was recently established in a rigorous way 
in Ref. 1351 . In this paper we generalize the results of the latter work to the case of an arbi- 
trary magnetic field. Using the equivalence between spin chains of HS type and inhomogeneous 
vertex models, in the case of spin 1/2 we are able to compute the large N limit of the chains' 
canonical partition function by means of the transfer matrix method. We obtain in this way an 
explicit unified expression of the free energy per site, which we systematically use for a de- 
tailed analysis of the thermodynamic properties. In particular, we deduce asymptotic formulas 
for the relevant thermodynamic quantities in both the zero magnetic field and the zero tempera- 
ture limits. Finally, we use the explicit formula for the free energy to show that spin chains of 
HS type are thermodynamically equivalent to a suitably weighted average of one-dimensional 
(homogeneous) Ising models. 
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As remarked throughout the preceding sections, some of the formulas for the thermodynamic 
quantities of the HS-type chains derived in this paper have previously appeared in the literature, 
especially in the case of the original Haldane-Shastry chain. It should be stressed, however, that 
our approach to the thermodynamics of these chains differs from previous related work in two 
essential ways. First of all, we treat all three families of spin chains of HS type in a unified way 
via their dispersion relation s(x) (cf. ( |4.17| i). As discussed in the previous section, this is possible 
due to the equivalence of these models (in the thermodynamic limit) to an inhomogeneous Ising 
model with couplings proportional to the values of the dispersion relation at the points x, = i/N, 
1 i : ^ N - 1 . On a more technical level, our method for studying the thermodynamics of spin 
chains of HS type is simply based on the evaluation of the canonical partition function in the 
thermodynamic limit through the standard transfer matrix approach. We thus bypass the more 
complex analysis relying on the equivalence of these models to an ideal gas of spinons used 
by Haldane to compute the thermodynamic functions of the original HS chain 11231 l43l . which 
requires the use of the grand canonical ensemble. 

The results presented in this paper can be generalized in two different directions. On the 
one hand, it should be possible to apply our analysis to general vertex models of the form ( |4.2| i, 
in which the coupling e, between consecutive spins is an arbitrary polynomial (or even a C 1 
function) in the variable x,, 1 ^ i ^ N — 1. These models, which have been recently studied 
by Basu-Mallick and collaborators ll35l |49l , share certain basic properties with spin chains of 
HS-type; for instance, their energies are normally distributed in the thermodynamic limit. On 
the other hand, it should also be of interest to extend our results to other types of HS chains, 
particularly those associated with the BCn and root systems (see, e.g., jl8 UTTI ). The first step 
in this direction would be to ascertain whether these spin chains are isospectral to suitable vertex 
models analogous to ( |4.2| i. Work on these and related topics, which is presently going on, shall 
be presented in subsequent publications. 
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Appendix A. Error term in the asymptotic formulas ( |5.22| i-( |5T23| 
Let us start with the error term in Eq. ( |5.22[ ). We must show that 



where we have set R - A/4. Note first of all that the LHS in the previous equation can be written 

as 





where the last term is easily seen to be 0(R 2 e R ). Hence we need only show that 
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Using the identity 



df 



R 2R) 



and integrating by parts we obtain 



f ^ f(i--r 1/2 -i 

Jo (1+e-OUV 

= 2R 1 - J I 

Jo I V R 2R 



g(t)e-'dt + 0(R 3 e- R ), 



where 



g(t) 



t z - It 2t z s 



2 Q -t 



(l+e-') 2 (l+e-') J 



The elementary inequality 



r, s s 

1 - Vl - s- - ^ — 
2 2 



implies that 



C R t 2 e-' \, fx-i/2 I 1 r 



df sC - I f 2 g(f)e dt + Q(R 3 e R ) = OiR' 1 ) , 



since g(t) = 0(f 2 ). 

Consider next the error term in Eq. ( |5.23[ >. We must now show that 



100 f 2 e- f 
o (l+e- f ) 2 



dt 



(r+iM 



dt 



(1 



or equivalently that 



M.J+DA t 2 e - t ^ I 

J„ (l+e-') 2 l ~\ l + y^J 



-1/2 



df = 0(A _i ), 



since 



r 



f 2 e-' 



i(y+lM (! + e ) 

Our assertion follows immediately from the elementary inequality 



1 -(1 +sy l/2 < -, s>0, 
2 



which implies that 

0^ 



<y+l)i { 2 e -t 



(l+e-') 



n2 



1-1 



4f 

r 2 ! 



-1/2 



df ^ 



2 f 3 e-' 



on account of the convergent character of the integral 



Jo (l+e-0 2 



df. 
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Appendix B. Computation of the definite integrals ( |6.25[ > 
Consider the definite integral 

( 1 + Vl + 4e<" x 



I(x 



;a)= | log 
Jo 



df. 



with a 2 — I. Performing the change of variable 



z= Vl +4e«') 



and taking into account that 



Li 



(I- V5 



log' 



f 1 + V5 ^ ;r 2 



15 



(cf. Ref. EH), we readily obtain 

Vl+4e«*) 



al(x; a) — | 



(1-V5) 



log(l-z) log(l-z) 



dz 



1 



log' 



Cl + Vl + 4e Q 



1 -z 
1 - Vl + 4e a 



Hence 



h =7(oo;-l)= — . 



As to the second integral in Eq. ( |6.25| l, we first note that 



h = lim \I(x\\)- 



From eq. (|B.1[) and the dilogarithm identity [44] 



-Li 2 (l-z)=Li 2 1 



1\ 1 



Z 2 



log z 



it follows that 



and therefore 



I(x; 1) = log' 



1 + Vl + 4e-' 



h = Li 2 (l) - — + lim 

13 x->oo 



+ Li, 1 - 



log' 



1 + Vl + 4e- 
1 + Vl + 4e- v \ 



= — + lim 

10 X~*00 



log' 



1 + Vl +4e x 



T 
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where we have made use of the equality Li 2 (l) = £(2) = n 2 /6 (cf. Ref. [44]). Since 



log 



( 



1 + Vl +4e v 



) 



2 



the limit in Eq. ( |B.3| l vanishes, and we finally obtain 



7T 



.2 



h = 



10 



(B.4) 
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